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Abstract 

Background: Mycobacteria can be quickly and simply identified by PCR restriction-enzyme analysis (PRA), but 
misidentification can occur because of similarities in band sizes that are critical for discriminating among species. 
Capillary electrophoresis can provide computer-aided band discrimination. The aim of this research was to develop 
an algorithm for identifying mycobacteria by combined rpoB duplex PRA (DPRA) and hsp65 PRA with capillary 
electrophoresis. 

Results: Three hundred and seventy-six acid-fast bacillus smear-positive BACTEC cultures, including 200 
Mycobacterium tuberculosis complexes (MTC) and 176 non-tuberculous mycobacteria (NTM) were analyzed. With 
combined hsp65 and rpoB DPRA, the accuracy rate was 100% (200 isolates) for the MTC and 91.4% (161 isolates) for 
the NTM. Among the discordant results (8.6%) for the NTM, one isolate of Mycobacterial species and an isolate of M. 
flavescens were found as new sub-types in hsp65 PRA. 

Conclusions: This effective and novel identification algorithm using combined rpoB DPRA and hsp65 PRA with 
capillary electrophoresis can rapidly identify mycobacteria and find new sub-types in hsp65 PRA. In addition, it is 
complementary to 16 S rDNA sequencing. 

Keywords: rpoB duplex polymerase chain reaction, hsp65 restriction fragment length polymorphism analysis .Capillary 
electrophoresis 



Background 

Detection and identification of mycobacteria in clinical 
specimens is a key issue in the therapy of pulmonary 
diseases because misidentification can lead to inappro- 
priate treatment. Traditionally, mycobacterial species are 
identified based on their growth rate, presence or ab- 
sence of pigmentation, and using biochemical assays of 
the isolates recovered from specimens. The biochemical 
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assays are time-consuming and labor-intensive, usually 
taking 1 to 2 months to complete, and assays for non- 
tuberculous mycobacteria (NTM) species can have poor 
reproducibility and provide ambiguous results [1,2]. 

By contrast, molecular identification, notably PCR- 
restriction enzyme analysis (PRA), is rapid and simple. 
The hsp65 PRA method, developed by Telenti et al. in 
1993, is a popular DNA-based method for mycobacteria 
identification [3]. Using hsp65 PRA, Wong et al. [4] 
reported 100% sensitivity and specificity in identifying 
Mycobacterium tuberculosis complexes but only 74.5% 
sensitivity in identifying NTM species. This misidentifica- 
tion may occur because of similarities in band sizes that 
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are critical for species discrimination [3]. An additional 
contributing factor is a lack of knowledge of all existing 
PRA profiles, especially among species that are very het- 
erogeneous, such as M. gordonae, M. scrofulaceum, and 
M. terrae complexes. Recently, capillary electrophoresis 
(CE) with computer analysis [5-9] has provided more pre- 
cise band discrimination than analysis by the naked eye. 

Previously, we developed an algorithm for mycobacte- 
rial species identification from acid-fast bacillus (AFB) 
smear-positive BACTEC tubes by combining the rpoB 
duplex PRA (DPR A) [10] and the key phenotypic char- 
acters of mycobacteria recovered from the tubes [11]. 
Using rpoB DPRA, we differentiated Mycobacterium 
tuberculosis complexes (MTC) from NTM with 235 base 
pair (bp) and 136 bp PCR amplicons in AFB smear- 
positive BACTEC cultures. The 136 bp rpoB duplex 
PCR amplicon was further digested with Mspl and 
Haelll (rpoB DPRA) to divide the NTM species into 
eight easily distinguishable groups (A-H) as described 
by Kim et al. [10]. Using two phenotypic characters 
(growth rate and photoreactivity on pigment production) 
and two simple biochemical assays (nitrate reduction test 
and Tween 80 hydrolysis test) [11], the mycobacterial spe- 
cies were identified. However, the sub-culture and bio- 
chemical tests for this algorithm took three weeks. 

In the present study, we developed a rapid and effect- 
ive algorithm for identification of mycobacteria by com- 
bined rpoB DPRA and hsp65 PRA with CE. 

Results 

Mycobacteria identification 

There were 376 AFB smear-positive BACTEC culture 
tubes (positive BACTEC cultures), including 200 MTC 
and 176 NTM-containing BACTEC cultures. A further 



20 bacteria were MGIT positive but AFB culture smear 
negative, and these were classed as contaminated and 
excluded from subsequent evaluation. 

By rpoB duplex PCR, all of the 200 MTC-containing 
BACTEC cultures and the 176 NTM-containing BACTEC 
cultures showed 235-bp and 136-bp PCR amplicons 
specific for MTC and NTM, respectively. The species 
were identified according to the flow chart shown in 
Figure 1. 

Concordant results from rpoB DPRA and hsp65 PRA 

Combining rpoB DPRA and hsp65 PRA with computer- 
aided CE gave an accuracy rate of 100% (200/200) for 
MTC and 91.4% (161/176) for NTM (Table 1). 

Discordant results from rpoB DPRA and hsp65 PRA 

There were 15 isolates (8.6%) of NTM with discordant 
results with rpoB DPRA and hsp65 PRA (Table 2). The 
two isolates, Mycobacterial species (A group) and M. 
flavescens (A group) identified by 16 S rDNA sequen- 
cing represented new patterns not available in the hsp65 
PRA databases and might be new sub-types in hsp65 
PRA. For Mycobacterial species, 16 S rDNA sequen- 
cing did not confirm the identity of the isolate but 
conventional biochemical identification showed it was 
M. mucogenicum. 

Development of a species identification algorithm 

The results in Tables 1 and 2 and the mycobacterial identi- 
fication flow chart (Figure 1) were used to develop a spe- 
cies identification algorithm by combining rpoB duplex 
PCR [10] and hsp65 PRA [3] using the most common 
74 patterns of 40 species in Table 3. In this algorithm 
(Table 3), we added M. gordonae types 3 and 4 to the 



Positive signals of MGIT cultures in the BACTEC 960 system 
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Figure 1 An flow chart of the identification of Mco bacteria I species from clinical specimens by combined rpoB duplex PCR(DCR) and 
hsp65 PCR-Restriction Fragment Length Polymorpgism analysis(PRA). 
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Table 1 Comparison of hsp65 RFLP, rpoB RFLP patterns, 16 S rDNA sequences and conventional biochemical 
identification in 361 isolates with concordant results 


rpoB RFLP pattern 


hsp65 RFLP pattern 


16 S rDNA sequence 
identification 


Conventional biochemical 
identification 


BACTEC culture Concordance 
number rate 


T 


M. tuberculosis type 1 


M. tuberculosis 


M. tuberculosis 


200 100%(200/200) 




NTM 


NTM 


NTM 


161 91.4%(161/176) 


A 


M. obscessus typel 


M. abscessus 


M. abscessus 


29 


A 


M. obscessus type 2 


M. abscessus 


M. abscessus 


41 


A 


M. fortuitum type 1 


M. fortuitum 


M. fortuitum 


33 


A 


M. fortuitum type 2 


M. fortuitum 


M. fortuitum 


2 


A 


M. peregrinum type 1 


M. peregrinum 


M. fortuitum* 


5 


A 


M. peregrinum type 2 


M. peregrinum 


M. fortuitum* 


8 


A 


M. peregrinum type 3 


M. peregrinum 


M. fortuitum* 


1 


A 


M. chelonoe type 1 


M. chelonae 


M. chelonae 


1 


A 


M. mucogenicum type 1 


M. mucogenicum 


M. mucogenicum 


2 


A 


M. smegmotis type 1 


M. smegmatis 


M. smegmatis 


2 


B 


M. avium subsp. avium type 2 


M. avium subsp. avium 


M. avium complex** 


2 


D 


M. kansasii type 1 


M. kansasii 


M. kansasii 


6 


D 


M. kansasii type 2 


M. kansasii 


M. kansasii 


1 


D 


M. kansasii type 6 


M. kansasii 


M. kansasii 


1 


D 


M. fr/V/a/e type 1 


M. triviale 


M. triviale 


1 


F 


M. malmoense type 1 


M. malmoense 


M. malmoense 


2 


F 


M. szu/ga/ type 1 


M. szulgai 


M. szulgai 


1 


F 


M. interjectum type 1 


M. interjectum 


M. interjectum 


1 


G 


M. intracellular type 1 


M. intracellular 


M. avium complex** 


14 


G 


M. gordonae type 1 


M. gordonae 


M. gordonae 


6 


G 


M. gordonae type 2 


M. gordonae 


M. gordonae 


1 


G 


M. gordonae type 5 


M. gordonae 


M. gordonae 


1 


Total 








361 



* M. peregrinum was identified as M. fortuitum by a conventional biochemical method. 

**M. avium subsp. avium and M. intracellular were identified as M. avium complex by a conventional biochemical method. 



F group. The accuracy rate for NTM by combining the 
two methods could reach 96.6% (170/176). 

Discussion 

Mycobacterial species identification is generally achieved 
by standard culture and biochemical methods [12]. Bio- 
chemical methods are labor-intensive and time-consum- 
ing, and not all are reproducible [2,13]. Recently, PR A 
and DPRA have been developed for molecular identifica- 
tion of mycobacterial species using different regions of 
hsp65, 16 S rDNA, 16 S-23 S rDNA spacer, dnaj, and 
rpoB as an amplification target [3,14-17]. The most com- 
mon method is hsp65 PRA, and 74 patterns for 40 spe- 
cies are available in the PRASITE database (http://app. 
chuv.ch/prasite/ index.html). 

Previous studies [18,19] have reported that hsp65 PRA 
is faster and more accurate for species identification than 



conventional (phenotypic or biochemical) testing. This is 
because more incorrect and ambiguous results are 
obtained with conventional methods. The results in our 
study (Tables 1 and 2) also support this finding. Incorrect 
and ambiguous results are caused by phenotypic homo- 
geneity among different species and phenotypic variability 
within species [18]. With by hsp65 PRA, some sub-species, 
such as M. kansasii, can be identified and rapid-growing 
mycobacterium can be divided into M. abscessus and M. 
chelonae, M. fortuitum and M. smegmatis [20], whereas 
these identifications are difficult with conventional methods 
[21]. As found in our study (Tables 1 and 2), M. peregrinum 
was identified as M. fortuitum and M. avium subsp, avium 
and M. intracellular were both identified as M. avium 
complex by the conventional biochemical method. 

However, hsp65 PRA limitations have been reported in 
some articles [22,23]. Failure to identify or incorrect 
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Table 2 Fifteen isolates of NTM species with discordant results from rpoB RFLP, hsp65 RFLP patterns, 16 S rDNA 
sequence and conventional biochemical identification 

No rpoB RFLP pattern hsp65 RFLP pattern 16 S rDNA sequence Conventional biochemical identification 

1 A fotEII : 242.8*, 214.0, 0 M. flavescens M. flavescens 

Hae\\\: 130.9, 140, 90.4, 49.7, 41.5, 37.1 



2 


A 


BstB\ :456.3, 0, 0 
Hae\\\:]92.6, 90.4, 82.0 


Mycobacterial species 


M. mucogenicum 


3 


D 


M. scrofuloceum type 1 


M. scrofulaceum 


M. scrofulaceum 


4 


G 


M. simiae type 5 


M. simiae 


M. simiae 


5 


G 


M. simioe type 5 


M. simiae 


M. simiae 


6 


F 


M. intracellular type 3 


M. intracellular 


M. avium complex** 


7 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


8 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


9 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


10 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


11 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


12 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


13 


F 


M. gordonae type 3 


M. gordonae 


M. gordonae 


14 


F 


M. gordonae type 4 


M. gordonae 


M. gordonae 


15 


F 


M. gordonae type 4 


M. gordonae 


M. gordonae 



^fragment size by CE. 

** M. avium subsp. avium and M. intracellular were identified as M. avium complex by a conventional biochemical method. 



identification of the species may occur because of simi- 
larities in band sizes critical for discriminating species, 
including difficult to distinguish M. tuberculosis complex 
(M. tuberculosis and M. bovis) [22], and closely related 
sub-species such as M. avium or M. gordonae, because 
of sequence heterogeneity [22]. In addition, technical 
problems can also cause misinterpretation or incorrect 
identification [23]. Patterns in PR A profiles are complex 
and difficult to interpret with the naked eye, especially 
when more detailed sub-types are included [21]. 

This study combined rpoB DPRA and hsp65 PRA to test 
both reference strains and clinical respiratory isolates. The 
mycobacterial identification flow chart (Figure 1) can iden- 
tify species to the sub-species level, and final species identi- 
fication can be obtained instantly with concordant results 
from the two PRA. M. gordonae has a highly variable gene 
sequence with 10 sub-types in hsp65 PRA, and there are 
two groups (G and F) in rpoB DPRA. Most M. gordonae is 
in the G group, but M. gordonae types 3 and 4 by hsp65 
PRA are in the F group (Tables 1 and 2). 

In addition, there were different rpoB DPRA results 
(Table 2) for M. simaie type 5 (G group but not E group), 
M. scrofulaceum type 1 (D group but not H group), and 
M. intracellular type 3 (F group but not G group). The 
identities of all of these isolates were finally confirmed by 
16 S rDNA sequencing. Variable numbers of restriction 
sites for Haelll in these species may derive from genetic 
sequence mutation. However, these species are included in 



the species identification algorithm even though they are 
uncommon isolates. 

Using the mycobacteria identification flow chart (Figure 1) 
and algorithm (Table 3), M. avium-intracellulare complex 
(MAC) can be easily divided into M. avium spp. avium and 
M. intracellular by both rpoB DPRA and hsp6S PRA. By 
contrast, this was not possible with the conventional 
method. Using the results in Table 3, some NTM species 
with identical or similar hsp65 PRA can be clearly grouped 
by rpoB DPRA (Table 4). Ambiguous results from hsp65 
PRA alone are easier to interpret with combined rpoB 
DPRA and hsp65 PRA. However, M. intermedium type 1 
and M. intracellular type 3 with identical hsp65 PRA and 
rpoB DPRA (G group) could not be differentiated further 
by this species identification algorithm and required 16 S 
rDNA sequencing for confirmation. 

Although 16 S rDNA sequencing is the standard method 
for mycobacterium species identification, it cannot differen- 
tiate some closely related rapid-growing mycobacterium 
species [24] or slow-growing M. kansasii and M. gastri that 
had identical 16 S rDNA sequences, but these can be differ- 
entiated by hsp65 PRA and rpoB DPRA. There are some 
reports [6,25] of conflicting results from different methods 
for mycobacterial species identification, probably caused by 
a failure of one method to identify all test strains correctly. 
Combining methods for mycobacterial species identifica- 
tion can improve the accuracy rate, avoid ambiguous 
results, and save time. 
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Table 3 A species identification algorithm 

nAlirmArnnicm analircic 

puiyiiiurpiiiMii diiaiybib 


by combining rpoB duplex PCR and hsp65 PCR-restriction fragment length 


rpoB DPCR-RFLP 


hsp65 RFLP 




Final species identification 


Pattern Msp I Haelll 


BstEII 


Hae III 


A 136 136 


440/0/0 


160/ 90 / 60/0 


M. vaccae type 1 




440/0/0 


160 / 85 / 55 /0 


M. flavescens type 3 




440/0/0 


140 / 55 / 50 /0 


M. flavescens type 1 




440/0/0 


130 / 115 / 70 / 60 


M. aurum type 2 




320 / 130 /0 


200 / 70 / 60 / 55 


M. immunogenum type 2 




320 / 130 /0 


200 / 60 / 55 / 50 


M. chelonae type 1 




320 / 130 /0 


145 / 70 / 60 / 55 


M. immunogenum type 1 




320 / 130 /0 


140/ 65 / 60 /0 


M. mucogenicum type 1 




320 /115 /0 


185 / 145 / 0 /0 


M. follox type 1 




320 /115 /0 


170/140 / 0 /0 


M. neoaurum type 1 




320 /115 /0 


145 / 65 / 60 /0 


M. mucogenicum type 2 




320 /115 /0 


140/ 90/60/0 


M. mucogenicum type 3 




235 / 210 /0 


200 / 70 / 60 / 50 


M. abscessus type 2 




235 / 210 /0 


180 / 135 / 70 / 50 


M. thermoresistibile type 1 




235 / 210 /0 


145 / 140 / 100 / 50 


M. peregrinum type 1 




235 / 210 /0 


145 / 70 / 60 / 55 


M. abscessus type 1 




235 / 210 /0 


140 / 125 / 100 / 50 


M. peregrinum type 2 




235 / 210 /0 


140/ 125 / 60 / 50 


M. senegolense type 3 




235 / 210 /0 


140/ 80 / 60 / 50 


M. phlei type 1 




235 / 210 /0 


130/ 80 / 60 /0 


M. celotum type 1 




235 / 130 / 85 


1 75 / 80 / 0 / 0 


M. aurum type 1 




235 / 130 / 85 


145 / 140/ 100/60 


M. peregrinum type 3 




235 / 130 / 85 


145 / 125 / 60 /0 


M. smegmotis type 1 




235 / 130 / 85 


140/ 125 / 60 / 50 


M. senegolense type 2 




235 / 120 / 85 


180/ 140 / 50 /0 


M. senegolense type 4 




235 / 120 / 85 


145 / 120 / 60/ 55 


M. fortuitum type 1 




235 / 120 / 85 


140/ 125 / 60/ 50 


M. senegolense type 1 




235 / 120 / 85 


140/ 120 / 60/ 55 


M. fortuitum type 2 




235 / 120 / 85 


135 / 90 / 85 /0 


M. fortuitum type 3 


B 136 108,28 


320 /115 /0 


140/ 90/60/0 


M. chitoe type 1 




235 / 210 /0 


145 / 130 / 0 /0 


M. avium subsp. avium type 3 




235 / 210 /0 


130 / 105 / 60 /0 


M. avium subsp. avium type 2 




235 / 210 /0 


130 / 105 / 0 /0 


M. avium subsp. avium type 1 




235 / 210 /0 


130 / 105 / 0 /0 


M. avium subsp. poratuberculosis type 1 


C 136 76,60 


235 / 120 / 85 


160/ 105 / 60 /0 


M. xenopi type 1 


D 75,61 or 136 


440/0/0 


170 / 130 / 0 /0 


M. triviole type 1 


75,57,4 


320 /115 /0 


130 / 95 / 75 / 60 


M. kansosii type 5 




235 / 210 /0 


190/ 105 / 80 /0 


M. ulcerans type 2 




(235 / 210 /0 


145 / 130 / 95 /0 


M. scrofuloceum type 1 *) 




235 / 210 /0 


145 / 105 / 80 /45/20 


M. marinum type 1 




235 / 210 /0 


145 / 105 / 80 /0 


M. ulcerans type 1 




235 / 210 /0 


130 / 105 / 80 /0 


M. kansasii type 1 
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Table 3 A species identification algorithm by combining rpoS duplex PCR and hsp65 PCR-restriction fragment length 
polymorphism analysis (Continued) 



75,61 or 
75,57,4 



108,28 



75,61 or 
75,57,4 



76,60 



75,61 or 
75,57,4 



76,32,28 



235 / 130 / 85 


140/ 105 / 70 /0 


M. shimodei type 1 


235 / 120 / 85 


130 / 115 / 75 / 60 


M. kansasii type 4 


235 / 130 / 85 


130 / 105 / 70 /0 


M. kansasii type 6 


235 / 130 / 85 


130 / 105 / 0 /0 


M. kansasii type 2 


235 / 130 / 85 


130 / 95 / 70 /0 


M. kansasii type 3 


440/0/0 


145 / 130 / 0 /0 


M. simiae type 5 


320 /115 /0 


185 / 140 / 0 /0 


M. terrae type 2 


320 /115 /0 


180/ 130 / 0 /0 


M. terrae type 1 


320 /115 /0 


145 / 130 / 0 /0 


M. simiae type 4 


320 /115 /0 


140/ 90/60/0 


M. nonchromogenicum type 2 


320 /115 /0 


140/ 60 / 50 /0 


M. terrae type 3 


320 /115 /0 


125 / 105 / 0 /0 


M. genavense type 1 


235 / 210 /0 


185 / 130 / 0 /0 


M. simiae type 1 


235 / 210 /0 


185 / 130 / 0 /0 


M. genavense type 2 


235 / 210 /0 


155 / 140 / 0 /0 


M. simiae type 2 


235 / 210 /0 


145 / 130 / 0 /0 


M. simiae type 6 


235 / 210 /0 


140/115 / 70/0 


M. terrae type 4 


235 / 130 / 85 


145 / 130 / 0 /0 


M. simiae type 3 


235 / 130 / 85 


130 / 105 / 70 /0 


M. gastri type 1 


235 / 120 / 85 


145 / 60 / 55 /0 


M. nonchromogenicumtype 1 


440/0/0 


130 / 105 / 70 /0 


M. szulgai type 1 


(320 / 115 /0 


130 / 115 / 60 /0 


M. gordonae type 4*) 


240/210/0 


130/110/0 


M. interject urn 


(235 / 210 /0 


145 / 130 / 0 /0 


M. intracellular type 3*) 


235 / 210 /0 


115/ 105/0/0 


M. asiaticum type 1 


235 / 130 / 85 


130 / 105 / 80 /0 


M. celatum type 2 


235 / 120 / 100 


145 / 105 / 80/0 


M. malmoense type 1 


235 / 210 /0 


145 / 105 / 80/0 


M. malmoense type 2 


(235 / 120 / 100 


130/ 115/0/0 


M. gordonae type 3*) 


(440 / 0 / 0 


145 / 130 / 0 /0 


M. simiae type 5*) 


320 /115 /0 


130/ 110 / 70/ 60 


M. gordonae type 8 


320 /115 /0 


130 / 115 / 60 /0 


M. gordonae type 4 


235 / 210 /0 


145 / 130 / 0 /0 


M. intermedium type 1 


235 / 210 /0 


145 / 130 / 0 /0 


M. intracellular type 3 


235 / 210 /0 


140/ 105 / 80 /0 


M. intracellular type 2 


235 / 210 /0 


130/ 115/0/0 


M. gordonae type 5 


235 / 210 /0 


120/ 115/ 110/0 


M. intracellular type 4 


235 / 130 / 85 


140/ 120 / 95 /0 


M. gordonae type 6 


235 / 120 /100 


160/115 / 60 /0 


M. gordonae type 9 


235 / 120 /100 


155/ 110/0/0 


M. gordonae type 7 


235 / 120 /100 


145 / 130 / 60 /0 


M. intracellular type 1 


235 / 120 /100 


130/ 115/0/0 


M. gordonae type 3 


235 / 120 /100 


130 / 110 / 95 /0 


M. gordonae type 1 0 


235 / 120 / 85 


160/115 / 60 /0 


M. gordonae type 1 
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Table 3 A species identification algorithm by combining rpoS duplex PCR and hsp65 PCR-restriction fragment length 
polymorphism analysis (Continued) 

235 /120 / 85 215 / 1 10 / 0 / 0 M. gordonae type 2 

H 75,61 or 66,60,10 235 / 210 / 0 145 /130 / 95 / 0 M. scrofuloceum type 1 

75,57,4 320 / 1 30 / 0 1 60 / 1 1 0 / 0 / 0 M. haemophilum type 1 

T 235 / 1 20 / 85 1 50 / 1 30 / 70 / 0 M. tuberculosis type 1 

(235 bp) 

(*)species from Table 2 were included in this algorithm. 



Many CE-based studies [5-9] in PCR-RFLP analysis 
have investigated improving band size discrimination. In 
one study by Chang et al [7], high-resolution CE gave 
more precise estimates of DNA fragment sizes than ana- 
lysis by the naked eye, and CE could detect low molecu- 
lar weight fragments (down to 12 bp). In our study, 
restriction fragments <50 bp were not selected for ana- 
lysis to avoid confusion with primer and primer dimer 
bands. We found that more hsp65 fragment differences 
than rpoB fragment (data not shown) may explain the 
size differences with highly variable sequence for species 
identification but difficult interpretation in hsp65 PRA. 

Some sub-types of NTM species are relevant to clin- 
ical management, such as the M. kansasii and MAC. M. 
kansasii type 1 is the most common type associated with 
human disease [26-28] because of its high pathogenicity. 
However, M. kansasii types 3-7 are most often isolated 



Table 4 Species with identical or similar hsp65 PRA but 
different groups in rpoB DPRA 



rpoB 
Group 


Species (type) 


hsp65 RFLP 
BsfEII 


Hae III 


A 


M. 


mucogenicum type 3 


320 /115 /0 


140/ 90/60/0 


B 


M. 


chitae type 1 


320 /115 /0 


140/ 90/60/0 


A 


M. 


mucogenicum type 2 


320 /115 /0 


145 / 65 / 60 /0 


E 


M. 


terrae type 3 


320 /115 /0 


140/ 60 / 50 /0 


A 


M. 


fa 1 lax type 1 


320 /115 /0 


185 / 145 / 0 /0 


E 


M. 


terrae type 2 


320 /115 /0 


185 / 140 / 0 /0 


A 


M. 


peregrinum type 2 


235 / 210 /0 


140 /125 /1 00/50 


H 


M. 


scrofulaceum type 1 


235 / 210 /0 


145 / 130 / 95 /0 


D 


M. 


kansasii type 6 


235 / 130 / 85 


130 / 105 / 70 /0 


E 


M. 


gastri type 1 


235 / 130 / 85 


130 / 105 / 70 /0 


F 


M. 


celatum type 2 


235 / 130 / 85 


130 / 105 / 80 /0 


D 


M. 


kansasii type 1 


235 / 210 /0 


130 / 105 / 80 /0 


F 


M. 


malmoense type 2 


235 / 210 /0 


145 / 105 / 80/0 


E 


M. 


simiae type 6 


235 / 210 /0 


145 / 130 / 0 /0 


G 


M. 


intermedium type 1 


235 / 210 /0 


145 / 130 / 0 /0 


G 


M. 


intracellular type 3 


235 / 210 /0 


145 / 130 / 0 /0 


F 


M. 


interjectum 


240/ 210 /0 


130/110/0 


G 


M. 


gordonae type 5 


235 / 210 /0 


130/ 115/0/0 



from the environment and rarely from humans, and 
have no significant role in clinical management [26,27]. 
MAC can be divided into M. avium subsp. avium and 
M. intracellular e because drug sensitivity test and clin- 
ical outcomes are different between these two sub- 
types [29,30]. 

It is important to identify NTM to the sub-type level 
both for epidemiologic data and for differentiating poten- 
tially pathogenic sub-types [26,27]. Combined rpoB DPRA 
and hsp65 PRA with capillary electrophoresis provides 
precise species identification and overcomes problems 
associated with discrimination by hsp65 PRA band sizes. 
This combined method takes around 2-3 days to 
complete in the laboratory once clinical isolates have been 
received. However, the identification algorithm has 
some limitations. First, it could not discriminate M. 
intermedium type 1 from M. intracellular type 3, and 
second, not every hospital laboratory will be equipped 
with the appropriate equipment for this method. 

Conclusion 

In conclusion, the novel flow chart and identification al- 
gorithm obtained by combined rpoB DPRA and hsp65 
PRA with capillary electrophoresis can easily differenti- 
ate MTC from NTM and identify mycobacterial species 
to the sub-type level, which is helpful for clinical man- 
agement. The results are complementary to 16 S rDNA 
sequencing, and the effective algorithm provides rapid 
and accurate mycobacterial species identification. 

Methods 

Mycobacterial isolates 

Fourteen mycobacterial reference strains including one 
MTC and 13 NTM strains and 376 clinical respiratory 
specimens, including sputum, broncho-alveolar lavage, 
and aspirated secretion from endotracheal tubes, were col- 
lected from January to July 2007 from Taichung Veterans 
General Hospital (Taichung, Taiwan). The respiratory spe- 
cimens were digested by a A/-acetyl-L-cysteine-NaOH de- 
contamination procedure, centrifugal concentration, and 
sputum dissolving agents [31]. The processed specimens 
or the concentrated specimens were inoculated into 
MGIT culture tubes and incubated in the BACTEC 960 
instrument at 37°C until a positive signal appeared. 
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Positive BACTEC cultures were smeared on glass slides 
and Kinyoun staining was used to screen for AFB [31]. 
Mycobacteria in the positive BACTEC cultures were iso- 
lated and identified by conventional methods [12,13] . 

Identification of mycobacterial isolates by conventional 
methods 

All 376 isolates from cultures in the AFB smear-positive 
BACTEC tubes were streaked on Lowenstein-Jensen (LJ) 
slants and incubated at 37°C under 5% C0 2 . Colonies on 
the LJ slants were used for species identification by con- 
ventional culture and biochemical methods [12,13]. These 
methods included growth rates, photoreactivity for pig- 
ment production, morphology in microcolonies on LJ 
slants, and biochemical tests, including nitrate reduction, 
arylsulfatase, Tween 80 hydrolysis, urease, semiquantita- 
tive catalase, tolerance to 5% NaCl and niacin production. 

Genomic DNA extraction 

Mycobacterial DNA was extracted from positive BACTEC 
cultures using a DTB specimen processing kit (Becton 
Dickinson, Franklin Lakes, NJ) according to the manufac- 
turer s instructions [11]. 

rpoB DPCR and rpoB DPRA 

The rpoB DPCR was performed using genomic DNA as 
template and primer pairs Tbcl (5'-CGTACGGTCGGCG 
AGCTGATCCAA-3>TbcR5 (5'-CCACCAGTCGGCGC 
TTGTGGGTCAA-3') and M5 (5'-GGAGCGGATGACC 
ACCCAGGACGTC-3>RM3 (5'-CAGCGGGTT GTTC 
TGGTCCATGAAC-3') as described by Kim et al. [10]. A 
235 bp DNA PCR amplicon from MTC and a 136 bp 
DNA PCR amplicon from NTM were specifically ampli- 
fied [10], and these two amplification products were ana- 
lyzed by electrophoresis on a 2% agarose gel (Seakem LE 
agarose, Cambrex, East Rutherford, NJ). 

For rpoB DPRA, the 136-bp DNA PCR amplicon was 
further digested with Mspl and Haelll after DPRA, and 
analyzed by electrophoresis on a 3% agarose gel (NuSieve 
3:1 agarose, Cambrex) or CE (eGene). The rpoB restriction 
fragment length polymorphism (RFLP) patterns were com- 
pared to eight groups described by Kim et al. [10]. 
Eight NTM reference strains (M abscessus ATCC 19977, 
M. avium subsp. avium ATCC 25291, M. kansasii ATCC 
12479, M. terrae ATCC 15755, M. szulgai ATCC 29716, 
M. intracellular ATCC 13950, M. scrofulaceum ATCC 
19981, M. xenopi ATCC 19250) from each rpoB group 
(A-H) were subjected to rpoB DPRA by CE (eGene). 

hsp65 PCR and hsp65PRA 

The hsp65 PCR was performed using genomic DNA as 
template and primer Tbll(5'-ACC AAC GAT GGT GTG 
TCC-3') and Tbl2 (ff-CTT GTC GAA CCG CAT ACC 
CT-3') as described by Telenti et al. [3]. A 439-bp DNA 



hsp65 PCR amplicon was specifically amplified from the 
extracted DNA, and the amplification product was ana- 
lyzed by electrophoresis on a 2% agarose gel (Seakem LE 
agarose, Cambrex). 

For hsp65 PRA, the 439-bp DNA hsp65 PCR amplicon 
was further digested with BstEll and Haelll after complet- 
ing hsp65 PCR, and analyzed by electrophoresis on a 3% 
agarose gel (NuSieve 3:1 agarose, Cambrex) or by CE 
(eGene). The sizes of the restriction fragment by hsp65 
PRA were compared to those reported on the PRASITE 
database (http://app.chuv.ch/prasite/index.html). Thirteen 
ATCC NTM reference strains and one MTC reference 
strain were subjected to hsp65 PRA by CE (eGene). 

eGene CE 

The restriction fragment sizes were estimated by the 
naked eye or by eGene CE according to the manufac- 
turers instructions using size markers (100 bp to 3 kb) 
and alignment markers (15 bp to 3 kb). 

Mycobacterial identification flow chart 

The mycobacterial identification flow chart is shown in 
Figure 1. 

16 S rDNA sequencing 

The 16 S rDNA sequencing of mycobacterial DNA as 
the reference standard method for mycobacterial species 
identification was carried out using primer pair 8FPL 
(5AGTTTGATCCTGGCTCAG 3') and 1492 (5'GGTT 
ACCTTGTTACGACT T 3') as described by Turenne et al. 
[32]. The species were identified by comparing the 16 S 
rDNA sequences with similar sequences from GenBank. 
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